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Abstract. M87 is the first known radio galaxy to 
emit very high energy (VHE) gamma-rays. During a 
monitoring program of M87, a rapid flare in VHE 
gamma-rays was detected by the MAGIC telescope 
in early 2008 Q). The flux was found to be variable 
on a timescale as short as 1 day, reaching 15% of the 
Crab Nebula flux above 350 GeV. In contrast, the flux 
at lower energies (150 GeV to 350 GeV) is compatible 
with being constant. We present light curves and 
energy spectra, and argue that the observed day- 
scale flux variability favours the M 87 core as source 
of the gamma-ray emission rather than the bright 
know HST-1 in the jet of M87. 
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I. Introduction 

M 87 is a giant elliptical radio galaxy (RG) of Fanaroff 
Riley I type (FR I; [2]) in the Virgo Cluster at a distance 
of 16 Mpc |3). It is powered by a supermassive black 
hole (BH) of (3.2 ±0.9) x 10 9 M Q 0. The M 87 jet was 
the first-ever observed [5|, and due to the proximity of 
M87, its morphological substructures can be resolved 
and a unique view of its innermost regions is possible. 
The jet, originating from the RG core, extends to 20" 
(El, equivalent to a 2 kpc projected linear distance). Sev- 
eral compact regions (knots) along its axis are resolved 
in the radio, optical, and X-ray regimes. These knots 
have similar morphologies in all wave bands, although 
the X-ray knots appear to be tens of pc closer to the 
core than the optical and radio knots Q. The variable 
brightness of the knots may be due to several shock 
fronts in the jet, being responsible for particle accel- 
eration and nonthermal emission. Superluminal motion 
of the knots has been observed in the optical (H and 
radio J3 wave bands, constraining the viewing angle of 
the jet to < 43° ± 4°. The innermost resolved bright 
knot HST-1 is located at 0.82" (64 pc) from the core 
and is the most prominent feature of the jet. HST-1 
has shown many flares exceeding the luminosity of the 
M87 core emission. Its X-ray brightness has increased 
by a factor 150 from 2000 to 2005 flU. A correlation 
between radio, optical, and X-ray luminosity points 
to a common origin of the emission. The measured 
superluminal motion in HST-1 is higher than in other 
knots, suggesting a viewing angle of < 19° for this 



part of the jet. The core itself is variable, too, and also 
shows a correlation between the emission levels from 
radio frequencies through X-rays [11 j. M87 was not 
detected by EGRET. The first hint of very high energy 
(VHE; > 250 GeV) 7-ray emission was reported by 
Aharonian et al. ( Ifl2l . 2003), and later confirmed by 
Aharonian et al. (O, 2006) and Acciari et al. ( fl4l . 
2008). The emission is variable on a timescale of years. 
The reported ~2 day variability ifTSl narrows down the 
size of the emission region to be on the order of the 
light-crossing time of the central BH. With its expected 
low accretion rate, the M 87 core radiation is not strong 
enough to attenuate significantly TeV 7-ray s even at 
5 Schwarzschild radii (Rs) away from the BH lfl5l . 
All this implies a production region in the immediate 
vicinity of the M 87 core. During later observations, no 
significant flux variation was found [ 14 1. An X-ray-VHE 
7-ray correlation is expected in most emission models, 
but was not unambiguously found so far. Whereas Aha- 
ronian et al. [13 1 claim a hint of a correlation between 
the soft (0.3-10 keV) X-rays at HST-1 and the VHE 
7-rays, Acciari et al. ( lfT4l 2008) find a year-by-year 
correlation between the 2-10 keV X-ray flux of the 
M 87 core and the VHE 7-ray emission instead, but do 
not observe a correlation between the two energy bands 
on shorter timescales. The radio to X-ray emission of 
the jet is due to nonthermal synchrotron radiation of 
relativistic electrons in the jet. The observed knots and 
flares in M 87 point to a complicated morphology with 
several shock fronts producing these electrons. While the 
majority of the currently known extragalactic VHE 7-ray 
emitters are blazars, M 87 is assumed to be a blazar not 
aligned to our line of sight fl6l . If the observed VHE 
emission from M 87 is associated with the innermost part 
of its jet, then blazar emission models may hold. Re- 
cently the H.E.S.S. collaboration reported the discovery 
of VHE radiation also from CentaurusA ifTTl . another 
FR I radio galaxy similar to M 87 under various aspects. 
In blazars without prominent disk or broad-line features, 
the VHE emission is explained by inverse Compton 
processes involving the synchrotron photons and their 
parent electron population (synchrotron self-Compton 
models; e.g., fl~8)). Alternatively, in hadronic models, 
interactions of a highly relativistic jet outflow with am- 
bient matter [19|[20|, proton-induced cascades ETI . or 
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Fig. 1: \ALPHA\ distribution for the overall data sam- 
ple. The background (gray histogram) is estimated using 
three OFF regions arranged symmetrically to the ON- 
source region with respect to the camera center. A 7-ray 
excess with a significance of 9.9 standard deviations is 
obtained. 



synchrotron proton radiation [22] [23] may produce VHE 
7-rays. In such a scenario, M 87 might also account for 
parts of the observed ultrahigh-energy cosmic rays n24ll . 
It should be noted that for M 87 the location of the VHE 
emission is still uncertain. Specific emission models for 
high-energy processes close to the core I25ll26ll27ll ll28l . 
in the large-scale jet [ 29 ) 1 30 1, or in the vicinity of a BH 
El (31] have been developed. The MAGIC collabora- 
tion performed monitoring observations of M 87 starting 
from 2008 January, sharing the task with the VERITAS 
and H.E.S.S. experiments. Here we report on MAGIC 
results from a subset of these data, revealing a variability 
timescale of M 87 of 1 day. The energy spectrum and 
light curves are discussed. 

II. Observations and data analysis 

The MAGIC telescope is located on the Canary Is- 
land of La Palma (2200 m above sea level, 28°45'N, 
17°54'W). MAGIC is a stand-alone imaging air 
Cerenkov telescope (IACT) with a 17 m diameter tessel- 
lated reflector dish. MAGIC has a low-energy threshold 
of 50-60 GeV (trigger threshold at small zenith an- 
gles). The accessible energy range extends up to tens 
of TeV with a typical energy resolution of 20%-30%, 
depending on the zenith angle and energy (32]. The 
data set comprises observations from 2008 January 30 
to 2008 February 11. These were performed in the 
wobble mode l33ll for 26.7 hr. The zenith angle of the 
observations ranges from 16° to 35°. After removing 
runs with unusually low trigger rates, mostly caused 
by bad weather conditions, the effective observing time 
amounts to 22.8 hr. The data were analyzed using the 
MAGIC standard calibration and analysis 113211 . The 
analysis is based on image parameters [34] [35] and the 
random forest (RF; [36]) method, which are used to 
define the so-called hadronness of each event. The cut 
in hadronness for 7/hadron separation was optimized 
on a contemporaneous data set of the Crab Nebula. 
After this cut the distribution of the angle ALPHA, 
which is the angle between the main image axis and 
the line between center of gravity of the image and the 
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Fig. 2: The night-by-night light curve for M87 as 
measured from 2008 January 30 (MJD 54495) to 2008 
February 1 1 (MJD 54507). The upper panel shows the 
flux in the energy bin 150-350 GeV, being consistent 
with a constant emission. The lower panel shows the 
integral flux above 350 GeV; flux variations are apparent 
on variability timescales down to 1 day. The inlay of 
the lower panel shows the light curve above 350 GeV 
in a 40 min time binning for the night with the highest 
flux (2008 February 1). The vertical arrows represent 
flux upper limits (95% c.l.) for the nights with negative 
excesses. 



source position in the camera, is used to determine the 
signal in the ON-source region. Three background (OFF) 
sky regions are chosen symmetrically to the ON-source 
region with respect to the camera center. The final cut 
\ALPHA\< 5° (Fig. 1) was also optimized on the Crab 
Nebula data to determine the number of excess events 
and the significance of the signal. The energies of the 7- 
ray candidates were also estimated using the RF method. 
To derive a differential energy spectrum, we applied 
looser cuts than those in Figure Q] to retain a higher 
number of 7-ray candidates and to lower the effective 
analysis threshold down to 150 GeV. Looser cuts also 
reduce systematic uncertainties between data and Monte- 
Carlo events, which is important for the estimation of 
the effective collection areas. The derived spectrum was 
unfolded to correct for the effects of the limited energy 
resolution of the detector J37). 

Finally, the spectrum and the light curves were cor- 
rected for trigger inefficiencies due to higher discrimi- 
nator thresholds during partial moon light and twilight 
conditions [38 1. These corrections are on the order of 
0%-20%. The data were also analyzed with an indepen- 
dent analysis yielding, within statistical errors, the same 
results. 
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Fig. 3: The differential energy spectrum of M87 for the 
total data sample. The horizontal error bars represent 
width of the energy bins. The best-fit function, assuming 
a power law, is given by the solid curve. The Crab 
Nebula spectrum is given by the dashed curve for 
reference. 
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Fig. 4: Differential energy spectra of M 87 divided into 
high (filled circles) and low (open circles) states. See text 
for the details. The best-fit functions, assuming power 
laws, are given by the black solid and red dashed-dotted 
curves, respectively. 



III. Results 

The \ALPHA\ distribution after so-called tight cuts 
is shown in Figure Q] The applied cuts are SIZE > 450 
photoelectrons and hadronness h < 0.02. After the final 
\ALPHA\ cut (resulting in an overall cut efficiency of 
37%), the total signal of 241 events over 349 normalized 
background events corresponds to an excess with a 
significance of 9.9 a along equation (17) in Li & Ma 
l39l . The highest flux was observed on 2008 February 
1 at a significance of 8.0 a. In searching for time 
variability, the data set was divided into nine subsets, 
one per observing night. In Figure ?? we show both 
the light curve above the analysis threshold (150-350 
GeV) and in the energy range at which MAGIC has the 
highest sensitivity for a variability search (>350 GeV). 
The low-energy range shows no significant variability 
with a xl of 12.6/8 (probability of P = 0.13) for 
a constant fit. Instead, in the energy range above 350 
GeV clear variability is found. A fit by a constant has 
a xl of 47.8/8 corresponding to P = 1.1 • 10~ 7 . The 
correlation coefficient between the two energy bins is 



r = -0.25~ 



0.33 



(1-cr errors), suggesting that there is no 



significant correlation, but we note rather large error bars 
in the low energy flux bin. We also investigated a night- 
by-night variability. There are five pairs of observations 
on consecutive nights in the total data set. We calculated 



individual probabilities Si for these pairs to have the 
same flux level and the corresponding significances. 
We then computed a combined significance S com b ll40l : 
Scomb = (%2Si)/y/fi, with n = 5. The resulting 
<Scomb = 5.6 cr, which we interpret as a proof that the 
flux varies on timescales of 1 day or below. Note that the 
1 day variability is claimed from this combined analysis 
rather than from the 2008 February 1 flare alone. We 
find our statistics not sufficient enough to determine the 
flare shape. Given the number of the observed changes 
in the flux level, the data belong to a complex of two, 
if not three, subflares. We also looked for shorter time 
variability, but in none of the observation nights there 
is a significant flux variation in the two energy bands. 
A typical example in a 40 minute binning is shown 
in the inset in Figure ?? for 2008 February 1. The 
averaged differential energy spectrum of M 87 (Fig. 3) 
extends from ^100 GeV to ~10 TeV and can be well 
approximated by a power law: 

-2.30±0.11 
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The errors are statistical only. We estimate an 11% 
systematic uncertainty in the normalization and 0.20 
for the spectral index J32). The measured values are 
in good agreement with the H.E.S.S. (spectral index 
T = -2.2 ±0.15; Aharonian et al. 2006) and VERITAS 
(r = -2.31 ± 0.17; Acciari et al. 2008) results. The 
observed spectrum is not significantly affected by the 
evolving extragalactic background light (EBL) due to 
the proximity of M87 fl5l . To investigate a possible 
hardening of the spectrum with increasing absolute flux 
level, we divided the data sample into high and low state 
subsamples. The high sample comprises the two nights 
with the highest flux above 350 GeV (February 1 and 8), 
while the low state comprises the nights of lower flux 
data (January 30, February 2, 4, and 11). Both the high 
and low state spectra (Fig. 4) can be well described by 
a power law: 

= /o 



dF 
dE 



ITeV/ 



10 



-12 



cm 2 s TeV 



with / high = (4.81 ± 0.82), r hi s h = (-2.21 ± 0.18) 
and /^ ow = (2.58 ± 0.92), T low = (-2.60 ± 0.30) for 
the high and low states, respectively. There is a marginal 
hardening of the spectral index with the higher flux on 
the level of 1-2 standard deviations, depending on the 
way the significance is calculated. This hardening is not 
significant, which can be a consequence of the fact that 
the two flux levels (states) differ by less than a factor 
of 2. 

IV. Discussion 

M87 is the first non-blazar radio galaxy known to 
emit VHE 7-rays and one of the best-studied extragalac- 
tic black hole systems. To enable long-term studies and 
assess the variability timescales of M87, the H.E.S.S., 
VERITAS, and MAGIC collaborations established a 
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regular, shared monitoring of M 87 and agreed on mutual 
alerts in case of a significant detection. During the 
MAGIC observations, a strong signal of 8 a significance 
was found on 2008 February 1, triggering the other 
IACTs as well as Swift observations. For the first time, 
we assessed the energy spectrum below 250 GeV, where 
our observations can be well described by a power 
law that shows no hint of any flattening. Our analysis 
revealed a variable (significance: 5.6 a) night-to-night 
7-ray flux above 350 GeV, while no variability was 
found in the 150-350 GeV range. We confirm the 
E > 730 GeV short-time variability of M 87 reported by 
Aharonian et al [13|. The observed variability timescale 
is on the order of or even below one day, restricting 
the emission region to a size of R < Ate 6 = 2.6 x 
10 15 cm = 2.6 SRs- The Doppler factor 6 is only 
relevant for an emission region not expanding while 
traversing the jet. In case of an expanding-jet hypothesis, 
the initial radius of the expanding shell is given by 
R* < cAt . The emission can occur very close to 
the BH, provided that the ambient photon density is 
low enough as to allow the propagation of VHE 7-rays. 
Otherwise the emission region must be located farther 
away from the BH. In the latter case, the variability 
constraints can be met only if the emitting plasma does 
not substantially expand while traversing the jet, or if 
it moves with 6 > 10. There exists no lower limit on 
the size of HST-1, and thus the flux variability cannot 
dismiss HST-1 as possible origin of the TeV flux. During 
the MAGIC observations, however, HST-1 was at a 
historically low X-ray flux level, whereas at the same 
time the luminosity of the M 87 core was at a historical 
maximum (D. Harris 2008, private communication). This 
strongly supports the core of M87 as the VHE 7-ray 
emission region. Our data alone cannot put strong con- 
straints on VHE 7-ray emission models. The relatively 
hard VHE spectrum found for M 87 (IT « -2.3) is not 
unique among the extragalactic VHE 7-ray sources if 
one considers intrinsic spectra, i.e., EBL corrected. Also, 
we did not measure a high-energy spectral cutoff. The 
found marginal spectral hardening may be interpreted as 
a similarity to other blazars detected at VHE, where such 
hardening has often been observed. Our results show that 
a dense TeV monitoring, as exercised by ground-based 
IACTs, has revealed highly interesting rapid flares in 
M87. This fastest variability observed so far in TeV 
7-rays observed in M87 restricts the size of the 7- 
emission region to the order of RS of the central BH of 
M87 and suggests the core of M87 rather than HST- 
1 as the origin of the TeV 7-rays. Results from the 
entire monitoring campaign, comprising data from other 
IACTs, will appear in a separate paper. 
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